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We investigate several models described by real scalar fields, searching for topological defects. 
Some models are described by a single field, and support one or two topological sectors, and others 
are two-field models, which support several topological sectors. Almost all the defect structures 
that we find are stable and finite energy solutions of first-order differential equations that solve the 
corresponding equations of motion. In particular, for the double sine-Gordon model we show how 
to find small and large BPS solutions as deformations of the BPS solution of the 0'' model. And 
also, for most of the two field models we find the corresponding integrating factors, which lead to 
the complete set of BPS solutions, nicely unveiling how they bifurcate among the several topological 
sectors. 
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I. INTRODUCTION 

This work deals with kinks and domain walls in sine-Gordon like models in (1, 1) spacetime dimensions. They are 
extended classical solutions of topological profile, supposed to play important role in several different contexts as for 
,. instance in condensed matter 0,13 and in high energy physics 
^ ' Models described by real scalar fields in (f , f ) space-time dimensions are among the simplest systems that support 
f — ■ topological solutions. Usually, the topological solutions are classical static solutions of the equations of motion, 
^\ \ with topological behavior related to the asymptotic form of the field configurations. The topological profile can be 
made quantitative, with the inclusion of the topological current, = (l/2)e'^''9i,0. It gives the topological charge 
Qt = (l/2)(0(a; oo) — (f){x — oo)), which is not zero when the asymptotic value of the field differs in both 
the positive and negative directions. To ensure that the classical solutions have finite energy, one requires that 
the asymptotic behavior of the solutions is identified with minima of the potential that defines the system under 
consideration, so in general the potential has to include at least two distinct minima in order for the system to 
support topological solutions. 

We can investigate real scalar fields in (3, 1) space-time dimensions, and now the topological solutions are named 
domain walls. These domain walls are bidimensional structures that carry surface tension, which is identified with the 
energy of the classical solutions that spring in (1, 1) space-time dimensions. The domain wall structures are supposed 
to play a role in applications to several different contexts, ranging from the low energy scale of condensed matter 
0,0,0,0 up to the high energy scale required in the physics of elementary particles, fields and cosmology 

There are at least three classes of models that support kinks or domain walls, and we further explore such models 
' in the next Sec. In the first class of models one deals with a single real scalar field, and the topological solutions 
are structureless. Examples of this are the sine-Gordon and (f)^ models ;3]- In the second class of models one also 
deals with a single real scalar field, but now the systems comprise at least two distinct domain walls. An example 
of this is the double sine-Gordon model, which has been investigated for instance in Refs. 
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the third class of models we deal with s yste ms defined by two real scalar fields, where one finds domain walls that 
admit internal structure [l^ ITEl IT^ ll^ llSl [l9l l20f . and junctions of domain walls, which appear in models of two 
fields when the potential contains non-collinear minima, as recently investigated for instance in Refs. |2ll . [2^ 0, 
HE HE m, HO, 1^ nil ' We study new possibilities in Sec. IIIII and there we investigate periodic systems, 
which present several topological sectors, which may bifurcate into richer structures. For most of them, we find the 
integrating factors, which lead us to all the BPS solutions the models engender. 

There are other motivations to investigate domain walls in models of field theory, one of them being related to 
the fact that the low e nerg y world volume dynamics of branes in string and M theory may be described by standard 
models in field theory |34l ISSl [s^ . Besides, one knows that field theory models of scalar fields may also be used to 
investigate properties of quasi-linear polymeric chains, as for instance in the appHcations of Refs. [33, Isl, |33, llfl , to 
describe solitary waves in ferroelectric crystals, the presence of twistons in polyethylene, and solitons in Langmuir 
films. 

Domain walls have been observed in several different scenarios in condensed matter, for instance in ferroelectric 
crystals 0], in one-dimensional nonlinear lattices 41j, and more recently in higher spatial dimensions - see |42j| and 
references therein. The potentials that appear in the models of field theory that we investigate in the present work 
are also of interest to map systems described by the Ginzburg-Landau equation, since they may be used to explore 
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the presence of fronts and interfaces that directly contribute to pattern formation in reaction-diffusion and in other 
spatially extended, periodically forced systems [2. Isl lisL li^. I45ll46l | . 

For completeness, in Sec. IIVI we investigate other topological sectors of the periodic two field models, where the 
domain wall solutions are of the non-BPS type. And finally, we end the paper in Sec. where we include our 
comments and conclusions. 



II. SINGLE FIELD MODELS 



In this work we are interested in field theory models that describe real scalar fields and support topological solutions 
of the Bogomol'nyi-Prasad-Sommerfield (BPS) type ^J. In the case of a single real scalar field 4>, we consider 
the Lagrange density 
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Here V{4>) is the potential, which identifies the particular model under consideration. We write the potential in the 
form 
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where W — W{(j)) is a smooth function of the field (p, and — dW/d4). In a supersymmetric theory W is the 
superpotential, and this is the way we name W in this work. 
The equation of motion for <f> = (j){x, t) has the general form 
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and for static solutions we get 



dx"^ 



W^W^^ (4) 
It was recently shown in Refs. |49l |50j that this equation of motion is equivalent to the first order equations 
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if one is searching for solutions that obey the boundary conditions Xvca.x^-oa'i'ix) = and XvccLx^-oo(d(\)l dx") — 0, 
where 0i is one among the several vacua {i^i, ^2, •■•} of the system. In this case the topological solutions are BPS (or 
anti-BPS) states. Their energies get minimized to the value V"^ = |AWij|, where AWij =Wi — Wj, with Wi standing 
for W{(f)i). The BPS solutions are defined by two vacuum states belonging to the set of minima that identify the 
several topological sectors of the model. 

In a recent work jSlj one has introduced a deformation prescription, in which one defines a deformation function 
/ — fi4')i which allows introducing new models presenting defect solutions. The deformation prescription goes as 
follows: if we change the potential V{(l)) to V {(/)), given by 

V{^)^V{ct>^ f)/{df/d^)' (6) 

where / = /(</>) is a well-defined, invcrtiblc function, then the defect solution of the new model can be obtained as 

0(x) = f-\cb{x)) (7) 

where (^(x) represents the defect solution of the former model. More recently, in another work one has enlarged 
the scope of the procedure, using other deformations. 
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A. Single field. One kink solutions 

We now turn attention to kinks and domain walls. Perhaps the most known example of this is given by the model, 
defined by the potential V{(j)) = (t;^ — 0^)^, where we are using natural units. In this model the domain wall can 
be represented by the solution (j)six) = ±ti tanh(Awx). The above potential can be written with the superpotential 
W{(j)) = A(w^0 — 0^/3), and the domain wall is of the BPS type. The wall tension corresponding to the BPS kink is 

Another example of system described by one field can be generated from the above 4>^ model, using the deformation 
method [sills^ with the deformation function f{(/)) — v sin(w^). In this case we obtain the sine-Gordon model defined 
by the potential 

V ^^X^cos^ivcj)) (8) 

with the kink and antikink solutions 

(f>{x) = ±— arcsin[tanh(Aux)] + kir (9) 

This potential can be written with the superpotential W{(l)) — (A/u) sin(i'0), and the kink is of the BPS type. There 
is an infinite set of minima given by (f> = kn/v, were k is positive or negative integer. The wall tension corresponding 
to the BPS wall is ts = 2|A/z;|. 



B. Single field. Two kink solutions 

We can build another class of models, in which the domain walls engender other features. We refer to models 
described by a single field, but the systems may now support two or more distinct topological solutions. An interesting 
example of this is the double sine-Gordon model defined by the potential, 

V = [cos\v(b) - a" sm\v(b)] ' (10) 

which can be generated from the cj)'^ model above, employing the deformation function fi{cf>) = ata.n{v(f>) or f2{<j>) — 
i tan(u0 — §)■ It is interesting to see that if one uses the deformation prescription with the /i-deformation, we get 
the solutions 

1 / 1 1 / X \A 

(pi[x) — ±— arctan — tanh(Ax) H (11) 

V \a J V 

However, if one uses the /2-deformation there are new solutions, given by 

, , , 1 , , / , NX (2n -|- 1)tt 

92[x) = =F- arctan (a tanh(Aa;)) H (12) 

V 2v 

These solutions represent large kinks and small kinks, which appear in the double sine-Gordon model. The novelty 
here is that we are finding the large and small kink solutions, using the deformation procedure of Refs. .SI. 52]. In 
Fig. 1, we depict the behavior of the potential in Ea. HlU|l . in terms of the parameter a. There we see the double 
sine-Gordon models, and we also show the topological sectors corresponding to the solutions (t>i{x) and (f>2{x). 

The potential (|10|l does not give double sine-Gordon behavior as A and v vary, for a = 1. Thus, since we shall 
concentrate the investigation on the double sine-Gordon model, which profile depends only on a, we prefer to change 
the potential to the form 

Vr{(j)) = [4r cos(0) 4- cos(2(^)] (13) 

where r is a parameter, real and positive. This potential is periodic, with period 27r, and for simplicity in the following 
we consider the interval ~2tt < <j) < 2tt. The value r = 1 distinguishes two regions, the region r G (0, 1) where the 
potential contains four minima, and the region r > 1, where the potential contains two minima. For r S (0, 1) the 
system supports two distinct wall configurations, the large wall and the small wall, which distinguish the two different 
barrier the model comprises in this case. The limits r — > and r — > 00 lead us back to the sine-Gordon model, with 
period tt and 27r, respectively. The double sine-Gordon model has been considered in several distinct applications, as 



FIG. 1: The double sine-Gordon model. The left figure shows the case X = v — 1, for a = 1 (dashed line), a = 1.5 (dash-dotted 
fine), and a = 2 (solid line). The right figure shows the case X = v = 1 for a = 1 (dotted line), a — 0.75 (dash-dotted line), 
and Q = 0.50 (solid line). 



for instance in Ref. [3.l9LllO|. where one investigates magnetic solitons in superfluid '^He, kink propagation in a model 
for poling in polyvinylidine fluoride, and properties related to the two different kinks that appear in such polymeric 
chain. More recently, it has also been used to model magnetic solitons in uniaxial antiferromagnetic systems |5,1 | . 
To expose new features of the double sine-Gordon model we rewrite Eq. in the form 

K(0) = -^[cos(0)-Kr]2 (14) 
1 + r 

where we have omitted an unimportant r-dependent constant. This potential is a particular case of the potential 
(|10|l . with a = — r)/(l -f r), v ~ 1/2, and A = \/l — r. Thus, it is a deformation of the model, and can be 
described by the superpotential 

l^((^) = -i=[sin(0) + r0] (15) 
Vl + r 




FIG. 2: The double sine-Gordon potential, depicted for r = 1/3 (dash-dotted line), r = 2/3 (dotted line), and r = 1 (solid 
line), to illustrate how the behavior of the model changes with r. Our notation is such that both V{(j>) and (j> are dimensionless, 
as we explain in the work. 

For r in the interval r G (0, 1) the minima of the potential are the singular points of the superpotential, dW/d(j) = 0. 
They are periodic, and for — 27r < (p < 2tt there are four minima, at the points (f) = ±7r± a{r), where a{r) = arccos(?'). 
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For r > 1 the minima are at ^ = ±7r, in the interval —2tt < 4> < 2tt. A closer inspection shows that for < r < 1 the 
local maxima at ±7r and the minima ±7r ± a{r) degenerate to the minima ±7r for r = 1, and remain there for r > 1. 
Thus, the parameter r induces a transition in the behavior of the double sine-Gordon model. 

The value r = 1 is the critical value, since it is the point where the system changes behavior: for r e (0, 1) this model 
supports minima that do not appear for r > 1. We illustrate the double sine-Gordon model in Fig. [2], where we depict 
the potential ^ for r = 1/3, 2/3 and for r = 1. Thus, it is very interesting to see that this model maps the presence 
of magnetic solitons in uniaxial systems |l53l |: in fact, the model is able to model both the antiferromagnetic, canted 
and weak ferromagnetic phases. In particular, in the canted phase there are two different domain walls connecting 
minima between two different barriers, as we show below. 

To get a better view of the the double sine-Gordon model, we examine the order parameter (j)(r), which is given 
by ±TT ± a(r) for < r < 1, so it goes continuously to ±n for r > 1. Also, the (squared) mass of the field can be 
obtained via the relation 



(16) 



where cj) is the corresponding minimum of the potential. For < r < 1 we get m^(r) — A — 4r, and for r > 1 we 
have m^{r) = 4(r — l)/(r -I- 1). We see that m{r) vanishes in the limit r ^ 1. These results indicate that r mimics 
a second order phase transition, a transition where the system goes from the case of two distinct phases to another 
one, engendering a single phase, as one illustrates in Fig. 3. 




FIG. 3: Plots of m^(r) (left) and 0(r) (right) for the double sine-Gordon potential, which illustrate how the behavior of the 
model changes with r. The quantities 0, m(r), and r that appear in these figures are all dimensionless. 



Let us first consider the case < r < 1. The energies of the BPS solutions are given as follows. For solutions 
connecting the minima — tt -|- a{r) and tt — a(r) the defect is large since it joins minima separated by a higher and 
wider barrier. We have 



TT — a{r) 



VT+7 

In the case of the minima tt — a{r) and tt -|- a{r) the defect is small and we get 

a{r) 



dsG 



4r- 



(17) 



(18) 



We notice that v^^q — t^^^ -I- ATTr/\/l + r, and the limit r — > 1 sends v^^q 



2^/2n and t^^Q 0, as expected. For 
the BPS states, from Eqs. Hll|l and (|12|l we can write the solutions explicitly. For instance, for solutions that connect 
the minima — tt -|- a{r) and tt — a{r) we get large kink solutions, which arc of the form 



i{x) — ±2 arctan 



1 1 ^ f 



(19) 
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For solutions that connect the minima tt ± a{r) and the minima — tt ± a(r) we get small kink solutions. They are 
given by 



4>s{x) — ±7r — 2 arctan 



(20) 



The case r > 1 is different. The minima are now at ±7r, and the model is similar to the standard sine-Gordon model. 
In this case we modify the potential of Eq. to 

VU) = — ^(cos^ ((/.) + 2rcos((/))+2r-l) (21) 
r + 1 

to make y(±7r) = 0. Thus, we can write this potential in terms of another superpotential, such that 

2 



= —^=^Jcos^{4>) + 2rcos((/)) + 2r - 1 (22) 
\/r + 1 

The presence of the square root complicates the calculation, and we have been unable to find explicit analytical 
solutions in this case. 

The potential in Eq. 114|l in the limit 7-^0 goes to Vb(0) = 1 + cos(20) which leads us back to the sine-Gordon 
model. Thus, we can suppose r small and use Vr{4>) to explore the double sine-Gordon model as a model controlled 
by a small parameter, in the vicinity of the sine-Gordon model. This feature is of direct interest to investigations 
that follow the lines of Ref. 54], which explores the vicinity of the BPS bound, and also in the case concerning the 
presence of internal modes of solitary waves, which seems to appear when one slightly modifies integrable models - 
see for instance Ref. |55l |. 

III. TWO FIELD MODELS 

We now turn attention to another class of models, which is described by two real scalar fields. In this case the 
domain walls may engender internal structure; see, e.g., Refs. jl4l Il5l IT6l| . In the present work we are interested in 
models of the sine-Gordon type, but now described by two real scalar fields. 

In the case of two real scalar fields (j) and Xi the potential is written in terms of the superpotential, in a way such 
that 



The equations of motion for static fields are 



V{<ly,x) = \wl + \wl (23) 



= W^W^^ + W^W^^ (24) 



ax 



W^W^^ + W^W^^ (25) 



which are solved by the first order equations 



Solutions to these first order equations constitute the BPS states of the model. They solve the equations of motion, and 
have energy minimized to t^^ = |AM^*-'| as in the case of a single field; here, however, AW'^^ = W{(j)i,Xi) ~^{4'jiXj)i 
since now we need the pair {(j)i,Xi) to represent each one of the vacuum states in the system of two fields. An 
important issue is that in the plane (0, x) one may have minima that are noncolinear, and this opens the possibility 
for junctions of defects. Moreover, in the case of two real scalar fields we can find a family of first-order equations 
that are equivalent to the pair of second-order equations of motion, but this requires that -I- W^x ^ the case 
of harmonic superpotentials [H^l • 

Models described by two real scalar fields are more intricate, and require more involved investigations. For instance, 
to solve the two first-order equations one can search for the integrating factor [5^ . However, the integrating factor is 
in general hard to find, and so we can use the trial orbit method, as suggested in [5^ in detail. 
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The models that we investigate are of interest both to ferromagnetic '53!| and ferroelectric systems. They are 
defined by the following general potential 



y{<P,x) ^v'^ {[r + cos{(f>) cos(p x)] ^ + / sin^ {<j)) sin^ {p x) } 



(27) 



where the parameter p is real and positive, with p 7^ 1 [the case p — 1 leads to single field models with a 7r/2 rotation 
in the (0, x) plane]. We will consider three distinct models, with r ~ 0, < r < 1, and r > 1. The case r = and 
r > 1 describe two coupled sine-Gordon models, and the case < r < 1 describe two coupled double sine-Gordon 
models. These models are new, and some of their features will be examined below. 
In general, all the three models can be described by the following superpotential 



W{(f>,x) ='n[r(l) + sin((/)) cos(px)] 



(28) 



and this will help simplifying the investigation, since the BPS solutions satisfy first order differential equations, which 
are simpler to solve compared to the equations of motion. 

In this section we search for BPS solutions in the models introduced in this paper. We investigate the first order 
equations corresponding to each one of the three models separately. For model 1 and model 2 we obtain the integrating 
factors, which allow solving the first order equations completely, finding all the family of orbits which connect the 
minima of the potential with energy minimized to BPS states. 

A. Model 1 

The first model is given by the potential H27() with r = 0. In terms of the superpotential (|28|l the first order equations 
are 



— ^ = 77 cos((/)) cos(px) 
ax 

^X 

— = -TIP sin((/))sin(px) 
ax 



(29) 
(30) 



There are minima at the points — ((2n+l)7r/2, m7r/p) a.nd Wn,m — {nn , (2to + 1) 7r/2p). In Fig. 4 we display 
these minima in the plane {4>,px)- 



px 



-2« 



FIG. 4: The minima of model 1 potential. The minima v axe represented by balls, and w by losangles. 

To obtain the most general solution to the first-order Eqs. (|29|l and (|30fl . we integrate the ordinary differential 
equation 



1 



cot(0) cot(px) 



(31) 
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which admits the integrating factor sin '^{px)i with <t = 1 + l/p^, and this determines all the orbits connecting the 
minima via the family of curves 



cos(0) = cisin-i(px) 



(32) 



where C is a real integration constant. In fig. 5, the behavior of a particular orbit with C, in the family H32() is 
displayed, where the critical value is C = p. 




FIG. 5: Orbits from model 1 given for: < C < p (dash-dotted line), C — p (solid line), and p < C < oo (dashed line). 



There are three type of BPS sectors connecting pairs of adjacent minima, namely: vvi for < C < p, vw for C — p, 
and VV2 for p < C < oo. The tensions of the BPS states are = tvv2 = 2|?7|, and = \rj\. Using the curve ^A2\ . 
the BPS equations (|29|l and H30|) decouple and can be rewritten as, for C G (0, oo) 



^ = r,cos(</))[l-i^2. ^Qg2. 
ax '■ 

^ = -r,psm{px)\l- E'^^" sin^^^ipx)]'' 
ax 



(33) 
(34) 



where D = 1/Cp, E — /p^ ^ and v = p^ . In general, we cannot obtain the analytical solutions of the equations l|33|) 
and H34|l . For this reason, we present some numerical solutions for each kind of orbit, in the following. 

Firstly, we investigate the existence of kinks in the vvi sector. For simplicity, we consider C = 0, which leads to 
orbits in the form of straight line segments connecting two vertical adjacent v minima, with (p = {2k + 1)tt/2. This 
reduces the potential H27|l to 

V ^^V^p^ sm\px) (35) 
The BPS equations (|29|l and (|30|) have the kink and antikink solutions: 

cf>^v,{x) = (2fc + l)| 

1 TT 

Xvviix) — ±- arccos(± tanh(?7p^a;)) + fc — 
P P 

In Fig. 6, we plot generic numerical solutions for < C < p. 

In the vw sector, the numerical solutions for C = 1/p are depicted in Fig. 7. 

Finally, we examine the vv2 sector. Two horizontal adjacent v minima are connected by the straight line orbit 
PX = kir, for C = oo. In this case, the potential (|27ll reduces to the form of Eq. Thus, the kink and antikink 
solutions of the BPS equations are 

Xvv2{x) = fc- (36) 



P 

<t>vv2 (2;) — ± arcsin(tanh(?7 x)) + kir 
In Fig. 8, we show numerical solutions in the more general case, for p < C < 00. 



(37) 
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FIG. 7: Numerical solutions for vw sector. 



Model 2 



Now we explore the second model, described by the potential lf?7|) with < r < 1. From the superpotential 
the first-order equations are 



^ = 77 [r + cos(0)cos(px)] 

— ^ = —r]psm{(j))sm{px) 
ax 



(38) 
(39) 



The set of minima are given by Vn.m = + (— l)"')7r/2 — Q:]/p + 2n7r, to vr/p) and Wn^m = {mr , ±[(1 + (— l)")7r/2 — 

a]/p + 2m7r/p), with a = arccos(r). In Fig. 9 we depict the minima in the plane ((j),px)- 

To obtain the most general solution to the first-order equations and we integrate the ordinary differential 
equation 



^ — — r- csc(</)) csc(px) cot((/)) cot(px) 

ax P P 



(40) 



which admits the integrating factor sin '^{px)i with a = 1 + 1/p'^ that determine all the orbits as the family of curves 



cos((/)) = F{x) = - sin'" \px) 
P 



C + r / sin "(px)rfx 



(41) 
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FIG. 8: Numerical solutions for VV2 sector. 



— ( 



FIG. 9: The minima of model 2 potential. The minima v are represented by balls and w by losangles. 

where C is a real integration constant. In Fig. 10, we depict behavior of a particular orbit with C, in the family l|41|l . 
where the critical values are Ci « —2.734 and C2 ~ 2.734, for p — 1/2 and r = 1/2. 

There are five different BPS sectors connecting pairs of adjacent minima. They are vvi for —00 < C < Ci , vwi 
for C — Ci, vv^ for Ci < C < C2, VW2 for C = C2, and t>t;2 for C2 < C < 00. The tensions of the BPS states are 



tyy-^ = 2\ri[ar — sin(Q;) 



I J tvV2 



2 |?7 [(tt ~ a)r + sin(a)]|, U 



oi^^vvi ~t~^ijD2)i tywi — 2 ^vv\: and tyy 



Using the curve (@lj, the BPS equation decouples, and can be rewritten in the form, for C ^ ±00, 

^ = -77Psin(px)[l-J^'(x)]* 



2 ^t;D2 • 



(42) 



We have been unable to obtain the general solutions of the above equation H42|) . For this reason, we present some 
numerical solutions for the several BPS sectors. 

Firstly, we investigate the existence of kinks in the vvi and vv2 sectors. Two horizontal adjacent minima are 
connected by the orbit — ^""1 for ~ —00 and C = 00. In this case, the BPS equations H29|) and (|30|l reduce to 



dx 
dx 
dx 



f] [r zL cos((/))] 




(43) 
(44) 
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FIG. 10: Orbits given for; — cxd < C < Ci (dasli-dotted line), C = Ci (solid line between dash-dotted and dotted lines), 
Ci < C < C2 (dotted line), C — C2 (solid line between dotted and dashed lines) C2 < C < 00 (dashed line). 



The potential l|27|l becomes 



l/± = i772[r±cos(0)]' (45) 



These potentials can be obtained from Eq. H10() with A = ^\/l — r^, v — 1/2, and a = y or y jz^, for V+ or V-, 

respectively. The solutions are given from Eq. (|11|) and Eq. H12f) . For the potential V+ the kink and antikink solutions 
are: 

Xvvi (x) = Xvv2 {x) = 2n- (46) 
P 



and 



and 



(f)^^^{x) = (2n+ l)7r±2arctanf\/^^ — -tanh(^Vl -r2 x)] (47) 

V 1 + 2 



(/) (x) = 2n7r ± 2 arctanfW i-t^ tanh(-yr^-r2 x)] (48) 

'■ V 1 — 2 ■' 

For the potential the procedure is similar. We get 

TT 

XvvAx) = Xvv2{x) = i^n + l)- (49) 

P 

and 

= 2n7r ± 2 arctan[i /| — - ta.nh{^ y/]~r^ x)] (50) 

and 



4>vv2{x) — (2n + l)7r ± 2 arctan[y^Y^~ tanh(^-\/l~^r^a;)] (51) 

The numerical solutions are presented in figs. 11 and 12, for —00 < C < Ci and C2 < C < 00, respectively. 
Kinks in the vwi, vv^ and vw2 sectors are displayed numerically in figs. 13, 14, and 15, for C = Ci, Ci < C < C2, 
and C = C2, respectively. 
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FIG. 11: Numerical solutions for vvi sector. 
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FIG. 12: Numerical solutions for vvi sector. 



Model 3 



The third model is given by the potential 127|l with r > 1. There is a set of minima at the points Un^m — 
(2m: , (2m + 1) n/p) and = ((2?^ + 1) tt , 2rm:/p). In Fig. 16 we show these minima in the plane {(j),px)- 

Here, the first-order equations are H88(l and (|39(l . Thus, the orbits are described by (|41|l . There are two kinds of 
BPS sectors connecting pairs of adjacent minima in this case, namely uu or vv for adjacent horizontal minima, and 
uv for adjacent diagonal minima. The tensions of the BPS states are i„„ = t^y = \ri\2Trr and t^v = 5 tw 

In general, we could not obtain analytical solutions of Eq. (|42|l . For this reason, we present some numerical solutions 
in the two BPS sectors. 

Kinks in the uu or vv sectors are found for C —^ 00, connecting two horizontal adjacent minima, with the orbit as 
a straight line px — kn. This reduces the potential (|27|l to the form H45() . In this case, the BPS equations H38() and 
(pTi^ have the solutions 



X{x) = k- 



(52) 



and 



fcTT ± 2 arctan 



/r+(-l) 



(-1) 



.tanhgV^ 



Ix 



(53) 



13 











_^ 

/ l|)(x) 


J 




a/p 




7t-a 


/ 

/ 


/X(x) 

/ 

/ 



X 



FIG. 13: Numerical solutions for vwy sector. 
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FIG. 14: Numerical solutions for vvi, sector. 

Particularly, for r = 1 the solutions take the forms: 

x(a;) = 2n— , <l)(x) = ±2 arctan(77a;)+A;7r (54) 
p 

or 

X(x) = (2n+ 1)-, 0(a;) = ±2 arccot(77x)+fc7r (55) 
p 

Some numerical solutions in the vu sector, for C = 10 and p = \/3/3, are depicted in Fig. 17. 

IV. NON-BPS SOLUTIONS 

In this section we search for non-BPS solutions in the models introduced in the former Sec. IIIII To identify the 
non-BPS sectors, we investigate the equations of motion corresponding to each one of the three models separately. 
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FIG. 15: Numerical solutions for VW2 sector. 
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FIG. 16: The set of minima of the model 3. 



A. Model 1 



From the potential H27() . with r = 0, we obtain the equations of motion 



d^cj) _ 1 
cix2 ~ 2 



if sin(2 0)[p2_(i+p2) cos2(px)] 



and 



= 1,72 sin(2px)[p' - (1 cos2(</>)] 



(56) 



(57) 



There are straight-line orbits connecting two adjacent vertical and horizontal w minima. The vertical orbit reduces 
the equations of motion H56() and H57|l to 



and 



dx2 







_ 1 2„3 -2, 



dx"^ 2 



7]^ p-^ sin''(2px) 



(58) 



(59) 
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FIG. 17: Numerical solutions for UU2 sector. 



The potential H27|l reduces to (jSU, with the kink and antikink solutions 



Xvviix) = ±- arccos[± tanh(77p^a;)] ' 
p 

The horizontal orbit reduces the equations of motion (|56|) and 157|) to 

<P4> 1 9 9 . /„ ,x 



fcTT 







The potential H27|l becomes V = sm'^{(f>), with the kink and antikink solutions 



XvV2 (•^) 



(2fc + 1) 



2p 



± arccos[± tanh(?7 ) + kir 



(60) 
(61) 

(62) 
(63) 



It is interesting to notice that if one considers the plane rotator model to describe the open states in the molecule 
of deoxyribonucleic acid (DNA), in Ref. 591 the author has obtained the coupled sine-Gordon equations 



i^- D^/Co)^ = 7^sin(w) + — sin(-u)cos(-w) 



p ^^'^^'") + ^ sm(-t;) cos(-u) 



with u = if + (p' and v = tp — ip' , where D, Co, I, Q and d are physical constants, and ip and tp' describe the angular 
displacements of two plane base-rotators. On the other hand, the equations of motions and (|77|l can be rewritten 
in the form 



d^ 

fx 

dx^ 



1 2 

2^^ 



:1 



i (1 - p2) sin(2 0) - 1 (1 + p2) sin(2 cj)) cos(2p x) 
i(l -p2)sin(2px) - ^(1 +p')sin(2px) cos(2(^ 



(64) 
(65) 



We notice that this model can be used as an alternative to the DNA model; this will be further examined in another 
work. 
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B. Model 2 

From the potential l(?7|) . with < r < 1, we obtain the equations of motion 

^-^ = -ry^ sin(0) Tp^ cos((/>) sin^(px) — cos(gx)['' + cos(p0) cos(gx)]l (66) 
dx'^ 2 ■' 

= -r/2psin(px)[p^cos(px)sin^(0)-cos((?!))[r + cos(0)cos(px)]] (67) 

There are straight-Hne orbits, (j) = kir, connecting two adjacent vertical w minima. The orbit reduces the potential 
(|?7jl to V± — [r± cos(px)] • For the potential V+, the solutions are 

where (jivv^i 2) given by H47|l and (|48|l . respectively. For the potential V^_, the solutions are 

0«;«.(i,2)(a;) = (2A: + l)7r, X«)tO(i,2) (a^) = %™(i,2) (a;) (69) 
where (pw^^ 2) given by H50|l and (|51|l . respectively. 

C. Model 3 

From the potential (|27(l . with r > 1, we obtain the equations of motion H66(l and (|67|l . There are straight-line orbits, 
(j) = kir, connecting two adjacent vertical it or minima. The orbit reduces the potential H27|l to 

V^^il^lr+i-l)" cosipxf. (70) 

The solutions are 

<P{x) = fcTT, x{x) = -Hx), (71) 

p 

where (f> is given by eq. (|53|l . 

The non-BPS solutions may be unstable, and may decay into stable BPS solutions. The decay of non-BPS states 
is out of the scope of this work, but it can be done following for instance the lines of Ref. W^. 

V. COMMENTS AND CONCLUSIONS 

In the present work we have investigated several models described by one and by two real scalar fields. The main 
investigations concern the search for BPS states, that is, for topological solutions that solve first order differential 
equations. These solutions minimize the energy to the Bogomol'nyi bound, which is given solely in terms of the 
superpotential, and the asymptotic value of the corresponding field configurations. 

The search for topological solutions is done at the classical level, and we have payed special attention to models 
introduced in Sec. In] some of them further explored in Sec. IIIII These investigations have shown how to write the 
double sine-Gordon model in terms of a superpotential, and to deal with the large and small kinks as BPS states, as 
solutions to first order differential equations. We have also investigated another model, decribed by two real scalar 
fields, defined by the potential of Eq. 127|l . It is similar to the model first investigated in Ref. ^6Q||, although here we 
deal only with periodic interactions. 

The classical investigations that we have developed are of interest in applications to nonlinear science, as for 
instance in the line of the work presented in Refs. [s^ IsM Is^ l40l| . where one uses field theory models to mimic 
nonlinear interactions in polymeric chains. For instance, the model 1 may be used as an alternative to the 0'* model 
that is standardly considered to mimic the polyacethylene (PA) chain, where Peierls instability appears due to single 
and double bound alternation in carbon atoms along the chain. The standard scenario leads to the very nice picture 
in which the distance carbon-carbon is mapped to the </)^ model with spontaneous symmetry breaking. The need for 
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spontaneous symmetry breaking is to reproduce the two degenerate states, which describe single-double and double- 
single bond alternations in the trans or zig-zag PA chain. Of course, this picture can become more interesting if one 
adds fermions to the system, via the standard Yukawa coupling. The alternative that we propose is to mimic the PA 
chain with model 1, which is very much similar to the ip'^ model at the classical level. We hope to explore this and 
similar ideas in the near future. 

A direct motivation that follows from the present investigations concerns the inclusion of fermions, to see how the 
fermions change the scenario we have just obtained. Another motivation concerns generalization of the field theory 
models to the case of complex fields, which may give rise to models engendering the continuum C/(l) symmetry. 
In this new scenario the models admit the presence of vortices, global and local, depending on the gauging of the 
global symmetry that appears when one changes the real field to a complex one. Furthermore, the two-field model 
investigated in Sec. IIIII is of direct use to describe interactions between two EEC's, and the corresponding vector 
solution describes the interface between the interacting condensates, as a good alternative to the recently proposed 
model of Ref. |67j . because in our model the center of the interface is asymmetric, containing differei it q uantities of 
each one of the two condensates that form the interface. This is more realist than the case studied in [^, where the 
center of the interface contains equal portions of each one of the two condensates. All the models here introduced 
are also of interest to pattern formation in spatially extended, periodically forced systems, governed by the Ginzb urg - 
Landau equation. This is a new route, different from the one proposed in [s^l, and further investigated in |3Ml39l l40| . 
which deals with solitons in ferroelectric cristals, in polyethylene, and in Langmuir films. The Ginzburg-Landau 
equation is appropriate for investigating fronts and interfaces, and their contributions to pattern formation. These 
and other specific issues are presently under investigation, and we hope to report on them in the near future. 
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